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ABSTRACT 

We present five new pulsating subdwarf B (sdB) stars discovered by the Kepler spacecraft 
during the asteroseismology survey phase. We perform time-series analysis on the nearly 
continuous month-long Kepler datasets of these 5 objects; these datasets provide nearly alias- 
^Jq \ free time-series photometry at unprecedented precision. Following an iterative prewhitening 

■ process we derive the pulsational frequency spectra of these stars, separating out artefacts 
\Q ' of known instrumental origin. We find that these new pulsating sdB stars are multiperiodic 

long -period pulsators of the VI 093 Her type, with the number of periodicities ranging from 8 
(KIC 8302197) to 53 (KIC 1 1558725). The frequencies and amplitudes are typical of g-mode 
pulsators of this type. We do not find any evidence for binarity in the five stars from their 
observed pulsation frequencies. As these are g— mode pulsators, we briefly looked for period 
spacings for mode identification, and found average spacings about 260 s and 145 s. This 
may indicate / = 1 and 2 patterns. Some modes may show evidence of rotational splitting. 
These discoveries complete the list of compact pulsators found in the survey phase. Of the 13 
K> ' compact pulsators, only one star was identified as a short-period (p-mode) V361 Hya pulsator, 

5_j \ while all other new pulsators turned out to be V1093Her class objects. Among the latter 

d ■ objects, two of them seemed to be pure V1093Her while the others show additional low 

amplitude peaks in the p— mode frequency range, suggesting their hybrid nature. Authenticity 
of these peaks will be tested with longer runs currently under analysis. 

Key words: stars: subdwarfs, oscillations (including pulsations) - space vehicles: instru- 
ments. 
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1 INTRODUCTION 

Hot subdwarf B (hereafter: sdB) stars are horizontal branch stars 
that consist of a helium rich core surrounded by a thin hydrogen 
envelope. The helium core sustains nuclear burning of helium into 
carbon and oxygen, though in some cases the core may have ex- 
hausted helium at the center. Their progenitors are low-mass main 

* E-mail:asb@iastate.edu 



sequence stars that are likely to have undergone the core helium 
flash (M < 2Mq). Once hydrogen in the core is exhausted, their 
progenitors evolve towards the tip of the red giant branch, where 
they lose all but a tiny fraction of hydrogen during (or prior to) the 
helium flash. Those stars that retain less than about 0.001 Mq of 
hydrogen lie on the high T e g end of the horizontal branch. This 
small fraction of hydrogen has an influence on the future evolution 
of sdB stars. As hydrogen shell burning cannot be sustained by the 
thin hydrogen envelope, the objects will move directly to the white 
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Figure la. Photometric data for KIC 7668647. Top panel shows the entire 
span of the data (with outliers removed); the bottom panel zooms in on part 
of the data to show representative variations. 



Figure lb. Same as Fig.[Ta] but for KIC 10001893. 



dwarfs cooling track, instead of going through the asymptotic giant 
branch phase. Several mechanisms for producing this small resid- 
ual hydrogen envelope t hat involve single-st ar or binary evolution 
have been proposed; e.g. lD'Cruz et al.ldl996l) and lHan et al. I fcOOl 
120031) . The hot subdwarfs have representative effective tempera- 
tures and surface gravities near 30000 K and log g ~5.5, respec- 
tively. For a general review of their observed properties, see Hebet 
(2009). 



Kilke nny et alj 1 19971) discovered the first hot subdwarf that 



showed stellar oscillations (hereafter: sdBV). The discovery pro- 
vided the first opportunity to use asteroseismology to look inside 
these enigmatic stars, with many more found with similar periods in 
the following years. No w, around 50 objects of this class are known 
dOstensen et al1l2010bl) . At the same time, independent theoretical 
work suggested that th ese stars could undergo nonradial pulsations 
(Charpinet et al. 1996), providing a "jump-start" to exploiting their 
asteroseismic potential. The first pulsators showed relatively short- 
period pulsations (the short-period sdBV stars, or V361 Hya stars, 
frequencies higher than about 2000 /iHz); these are attributed to 
pressure modes (p-modes) and are driven in the outer part of the 
star. As discovered later, sdB stars can also s how variation wit h 
pulsation periods an order-of-magnitude longer dGreen et al.l2003l) . 
The variability in these long period sdBV stars, (or V1093Her 
stars) show frequencies below about 500/rfIz; this variability is 
caused by gravity modes (g -modes) as suggested by theoretical 
models (Fontaine et al. 2003|). The g-modes are established deeper 
in the stars than the pressure modes. Hybrid sta rs, pulsating in bot h 
kind of modes were also found observationally dSchuh et al.l2006t) . 
As they exhibit both types of modes, models of these stars can give 
better constraints on their interiors since the pulsations probe both 
the deeper layers as well as layers nearer the surface. 

Observations of long period sdBV stars, including detection 
of new pulsators, are not easy from the ground. Relatively long 
period photometric variations (with only a few cycles per night) 
can be affected by variable sky transparency. To make defini- 
tive determination of the pulsation frequency requires extended 
photometric campaigns, preferably at several sites widely spaced 
in longitude to reduce day/night aliasing. However, Kepler pro- 
vides long-duration continuous, homogenous, evenly spaced time- 
series photometry, making it an ideal instrument for asteroseismol- 



ogy. The Kepler science g oals, mission design, and overall per- 
for mance are reviewed b y lBorucki et al.l (12010). Koch et al. 
and lJenkins et al.l d20ld) . 

Objects classified as hot subdwarfs located in field of view 
of the Kepler photometer were observed during the survey phase 
in the first year of the mission. As the spacecraft is rolled every 
3 months, the survey phase was divided into 4 quarters. Analysis 
of the observations of compact stars from the first half of the sur- 
vey phase (Q0, Ql and Q2) is de scribed in the following papers: 
Osten sen et aUl2010d. Paper I). dKawaler et alJl20ld3 . Paper II), 



Reed et alj|20ld. Pap er III), dVan Grootel et alihoid Pap er IV) 
Kawaler et alj2010bl . Paper V) and lOstensen et allfeOlOcj) . 



In this paper we focus on pulsating sdB stars discovered 
with data from the last 6 months of the survey phase (Q3 and 
Q4). These are (ordered by quarters): KIC 7668647, KIC 8302197, 
KIC 10001893, KIC 10553698 and KIC 11558725. In the Q3 data 
we found 4 sdBV of V1093 Her class, while in the Q4 dataset we 
found a fifth object. All stars presented in this paper are spectro- 
scopically confirmed sdB stars (see Paper VI). In Table [2] we pro- 
vide their spectroscopic parameters. They were obtained by fitting 
spectra to model grids in order to derive effective temperature, sur- 
face gravity and photospheric helium abundance, which was an in- 
dication whether an object is a hot subdwarf or not. Spectroscopic 
data along with fitting procedure are p rovided in the leading paper 
for the final half of the survey phase (Ost ensen et alj|201 ll Paper 
VI). 

Although we did not find any predominantly p-mode pulsators 
from Q3 or Q4, we detected many peaks in the short period regime 
in many targets. However, these are few, compared with the large 
number of peaks identified in the longer period (/-mode range. 
Moreover, Kepler data have a variety of instrumental issues caus- 
ing some artefacts to appear in the photometric data at frequencies 
comparable to the p-modes seen in sdBV stars. Some of these arte- 
facts are related to the long-cadence (LC) readout timings, are well- 
characterised, and their frequencies are known to high precision (as 
described in details in the Kepler Data Release Notes and Gilliland 
et al. 2010). We believe, though, that some single peaks in the short 
period regime may also be associated with as-yet unidentified arte- 
facts, since they are seen in more than one object. Based on this, we 
treat all peaks with frequencies higher than 1000/iHz with caution. 
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Table 1. Log of observations. 
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Figure lc. Same as Fig. [la] but for KIC 8302197. Here, the variability is 
not as obvious. Compared to all stars in this sample, the average noise level 
is the largest, and amplitudes of detected peaks are the smallest. 

On the other hand if they are confirmed to be intrinsic to the stars, 
then all stars presented in this paper are hybrid pulsators. 

We have also found peaks in the "transition region" between 
the p- and g-modes (i.e. between 500 and 2000 /iHz). These peaks 
we re observed in the s dB V stars Balloon 090 1 00001, KIC 2697388 
by iBaran et al] j2009t) and [Reed et al.S d2010h . respectively. Their 
authenticity also needs confirmation. 

The overall survey strategy and properties of the candidates 
have been described in prior papers; Paper VI presents all selected 
targets for Q3 and Q4. Further analysis of the pulsation spectra 
of these st ars (and several fr om the first part of the survey) is de- 
scribed by (Reed et alj20Tl Paper VIII), who attempt to disentan- 
gle asymptotic behaviour of p-mode pulsations in V1093 Her class 
stars observer by Kepler and CoRoT. 



Figure Id. Same as Fig.^a] but for KIC 1 1558725. 
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2 OBSERVATIONS 

Data from Q3 cover 3 months, starting from BJD = 2455092.7 123 
(Q3.1) and ending at 2455182.0065 (Q3.3). Q4 covers the next 3 
months from BJD = 2455184.8679 (Q4.1) to 2455274.7137 (Q4.3). 
One safe-mode event occurred during Q3. As a result, data for two 
days (between BJD = 2455154 and 2455156) are not available. As 
it happened at the beginning of the third month in Q3, it did not 
cause any gaps in the data presented here. However, the third month 
of Q3 (Q3.3) was truncated by those 2 days, slightly reducing the 
frequency resolution as well as slightly increasing the noise level in 
the amplitude spectrum as a result of less data (and fewer pulsation 



Figure le. Same as Fig.[Ti] but for KIC 10553698. 

cycles). Two stars in this paper are affected by the safe-mode event: 
KIC 8302197 and KIC 1 1558725. 

Beside this safe-mode event, the Kepler spacecraft lost its 
fine pointing control twice. The first was a brief interval be- 
tween BJD = 245099.91 and 2455100.04, which was quickly cor- 
rected. The second, longer pointing control error fell between 
BJD = 2455 113.05 and 2455113.83, and required correction via 
ground command. Only data on KIC 7668647 were affected by this 
fault. A log of the photometric data is presented in Table Q] 



© 2010 RAS, MNRAS 000.ITHT51 



4 A. S. Baran et al. 



Table 2. Spectroscopic parameters with their errors of five stars analyzed 
in this paper. 



KIC 


K P 


T eS [K] 


logs 


log N He m H 


7668647 


15.4 


27,700(300) 


5.45(4) 


-2.5(1) 


10001893 


15.8 


26,700(300) 


5.30(4) 


-2.9(1) 


8302197 


16.4 


26,400(300) 


5.32(4) 


-2.7(1) 


11558725 


14.9 


27,400(200) 


5.37(3) 


-2.8(1) 


10553698 


15.1 


27,600(400) 


5.33(5) 


-2.9(2) 



All data were obtained thro ugh the Kepler Astero seismic Sci- 
ence Consortium (KASC) - see lGilliland etafl faOlCh . The KASC 
data include raw as well as calibrated fluxes. Calibration takes into 
account all known instrumental issues as well as contamination 
from nearby objects. Since spectroscopic parameters of stars in the 
Kepler database are, generally, not precise, we believe that this con- 
tamination may not be evaluated correctly. That is why we decided 
to use raw fluxes for all stars in our analysis and not to account 
for contamination. This may lead to decreased amplitudes of de- 
tected modes by a factor of (1 — c), where c is a contamination 
factor. Future processing, supported by ground based observation 
of the relevant fields, will allow more accurate determination of the 
correction factor. 

We analyzed short cadence (SC) data. One SC data point is a 
sum of nine integration of 6.02 s followed by 0.52 s readout result- 
ing in 58.85 s time resolution. This in turn gives 8496.356 /iHz as 
a Nyquist frequency. Preparing data for Fourier analysis, first we 
removed all SC points flagged with errant times or fluxes. Then, all 
remaining points (if necessary) were subject to de-trending. This 
was done via a cubic spline fit calculated on 0.3 days subsets and 
subtracted from original data. De-trending with this interval will 
suppress any frequencies below 23 /iHz. In the next step, we re- 
moved all outliers using 4a clipping. Only then, the data were sub- 
ject to Fourier analysis. All reduced data for each star, with shorter 
representative segments of the light curves, are presented in FigfTal 
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Figure 2a. The amplitude spectrum of KIC 7668647. The top panel shows 
theoriginal amplitude spectrum up to the Nyquist frequency, while the sec- 
ond panel from the top shows the low frequency domain. The third panel 
shows the residual amplitude spectrum after removal of the 9 peaks with 
the highest amplitudes (those frequencies are marked with short dashed 
lines). The bottom panel shows the Fourier transform of the residuals af- 
ter 20 peaks have been removed. Short solid lines denote the additional 1 1 
periodicities removed from the data. The horizontal line represents S/N = 4 
threshold at 0.068 mma. 



3 AMPLITUDE SPECTRA 

We used Fourier analysis to identify candidate pulsation frequen- 
cies. We then did a nonlinear least squares fit including each pe- 
riodicity in the form Ai sin(ujit + cj>i). This was done in an iter- 
ative process, following the standard pre-whitening procedure and 
continued until all peaks with amplitudes satisfying signal-to-noise 
ratio (hereafter: S/N) ^ 4 rule had been removed. Further details 
and examples can be found in Papers II, III, and V. 

KIC 7668647 The amplitude spectrum of this object is domi- 
nated by a number of peaks between 100 and 400 /iHz and is shown 
in Fig.|2a| The top panel presents all frequencies up to the Nyquist 
frequency, while the second panel isolates the frequency range be- 
tween 100 and 400 /iHz. Peaks located in this frequency range re- 
sult from g-mode pulsation. We can easily identify more than ten 
peaks in the original amplitude spectrum. The amplitudes less than, 
or on the order of, 1 mma are typical for V1093 Her stars. 

Aside from these peaks, the top panel of Fig.|2alshows several 
peaks at higher frequencies. These are signatures of the artefacts 
at multiples of the LC frequency of 566.4 pRz (7 th , 8 th and 14 th 
harmonics) and two other peaks at about 1000 and 4700 /iHz. The 
peak at 1000 /iHz is located in a "transition region" where peaks 
in other sdBV stars (mentioned in Section Q3 were identified, in- 



cluding analysis of ground based data for Balloon 090100001). The 
frequency of 4700 /iHz for the other peak is typical for p-modes in 
sdB pulsators. We should note however that such single peaks are 
observed for almost all sdBV stars detected with Kepler photome- 
try, and the ones here were barely detected in our data. Although 
they are not associated with any known artefacts, they should be 
treated with caution now. With more data coming (runs of several 
months up to one year) we will be able to say with greater certainty 
whether or not they are intrinsic to the star. If confirmed, this object 
will be exposed as a hybrid pulsator. 

These two peaks were taken into account in our analysis as 
they are isolated and do not affect determination of the parameters 
of the lower-frequency g-modes. In total we detected 20 periodic- 
ities with S/N^4. All pre-whitened frequencies are listed in Table 
[3] The bottom panel of Fig. [2a] still contains peaks with significant 
amplitude. However, all those peaks lie close to already-removed 
periodicities. This may indicate amplitude/phase variation of the 
periodicities, or that there are a number of nearly-degenerate oscil- 
lation modes that remain unresolved in this 30-day observation. We 
did not attempt to pre-whiten those residual peaks. The mean noise 
level in the final residual amplitude spectrum is 0.017 mma. 

KIC 10001893 As in the case of the previous object, the am- 
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Table 3. Results of the pre-whitening process for KIC7668647. The phases are given at mean epoch BJD = 2455108. 261682. The numbers in parentheses are 
the errors of the last digits. In all tables presenting list of detected peaks, horizontal lines divide whole range of frequencies into three regions: g-mode region 
(defined by the bulk of peaks), "transition region" between g-modes and 2000 /iHz and frequencies - above 2000 /^Hz. 



ID 


Frequency [/iHz] 


Period[s] 


Amplitude [mma] 


phase [rad] 


S/N 


fl 


109.963(39) 


9094.0(3.2) 


0.074(14) 


2.11(19) 


4.3 


h 


115.878(21) 


8629.8(1.6) 


0.138(14) 


4.59(10) 


8.1 


f 3 


119.514(10) 


8367.2(0.7) 


0.307(14) 


1.234(45) 


18.0 


t4 


131.477(31) 


7605.9(1.8) 


0.093(14) 


5.54(15) 


5.5 


fs 


145.709(3) 


6862.99(0.14) 


0.922(14) 


4.562(15) 


54.2 


f 6 


159.026(31) 


6288.3(1.2) 


0.092(14) 


6.01(15) 


5.4 


fr 


170.130(8) 


5877.86(0.28) 


0.347(14) 


5.774(40) 


20.4 


fo 

E 8 


1 89 076H4^ 


S988 QH frt 




o ^nn6~i 


5 1 


f 9 


194.522(5) 


5140.81(0.13) 


0.640(14) 


1.802(22) 


37.6 


fio 


196.084(10) 


5099.86(0.26) 


0.289(14) 


2.769(48) 


17.0 


fu 


201.319(36) 


4967.2(1.0) 


0.080(14) 


3.16(17) 


4.7 


fl2 


204.513(18) 


4889.66(0.43) 


0.163(14) 


4.70(9) 


9.6 


fl3 


226.532(8) 


4414.39(0.16) 


0.372(14) 


4.648(37) 


21.9 


fl4 


241.143(35) 


4146.9(0.6) 


0.083(14) 


4.12(17) 


4.9 


fl5 


243.157(24) 


4112.57(0.41) 


0.123(14) 


4.70(11) 


7.2 


fl6 


261.736(28) 


3820.64(0.41) 


0.103(14) 


2.90(14) 


6.1 


hr 


295.183(16) 


3387.73(0.18) 


0.185(14) 


2.66(8) 


10.9 


hs 


345.795(16) 


2891.89(0.13) 


0.187(14) 


2.81(8) 


11.0 


fl9 


1043.00(37) 


958.77(0.03) 


0.079(14) 


2.41(18) 


4.6 


f"2() 


4737.812(31) 


211.068(0.001) 


0.095(14) 


2.12(15) 


5.6 



Table 4. Results of the pre-whitening process for KIC 10001893. The phases are given at mean epoch BJD = 2455139. 431910. The numbers in parentheses 
are the errors of the last digits. 



ID Frequency[^tHz] Period[s] Amplitude [mma] phase [rad] S/N 



fl 


77.522(24) 


12899.5(4.0) 


0.157(18) 


4.44(12) 


7.2 


h 


114.370(32) 


8743.6(2.5) 


0.120(18) 


2.04(15) 


5.4 


h 


125.910(30) 


7942.2(1.9) 


0.128(18) 


2.82(14) 


5.8 


U 


130.551(31) 


7659.8(1.8) 


0.125(18) 


5.90(15) 


5.7 


h 


135.338(24) 


7388.9(1.3) 


0.159(18) 


3.99(11) 


7.2 


fe 


140.592(38) 


7112.8(1.9) 


0.101(18) 


1.51(18) 


4.6 


f'7 


146.087(13) 


6845.2(0.6) 


0.297(18) 


3.71(6) 


13.5 


f 8 


152.054(7) 


6576.61(0.30) 


0.593(18) 


1.716(31) 


27.0 


fg 


158.001(39) 


6329.1(1.6) 


0.100(18) 


5.49(18) 


4.5 


tin 


165.301(34) 


6049.6(1.2) 


0.113(18) 


4.21(16) 


5.1 


fll 


172.628(22) 


5792.8(0.8) 


0.179(18) 


1.94(10) 


8.1 


fl2 


180.499(32) 


5540.2(1.0) 


0.121(18) 


5.65(15) 


5.5 


fl3 


189.072(20) 


5289.0(0.6) 


0.197(18) 


1.30(9) 


8.9 


fl4 


199.112(7) 


5022.30(0.18) 


0.575(18) 


5.721(32) 


26.1 


fir, 


204.713(28) 


4884.9(0.7) 


0.136(18) 


5.76(13) 


6.2 


fie 


210.688(6) 


4746.35(0.14) 


0.644(18) 


4.233(28) 


29.3 


fl7 


255.153(20) 


3919.22(0.31) 


0.189(18) 


6.14(10) 


8.6 


fl8 


262.954(9) 


3802.95(0.13) 


0.426(18) 


4.039(42) 


19.4 


fig 


274.301(4) 


3645.63(0.05) 


1.023(18) 


3.512(18) 


46.5 


fzo 


286.025(23) 


3496.20(0.28) 


0.168(18) 


3.09(11) 


7.6 


hi 


298.656(16) 


3348.33(0.18) 


0.234(18) 


2.63(8) 


10.6 


f'22 


323.988(5) 


3086.53(0.05) 


0.782(18) 


4.452(23) 


35.6 


h3 


359.685(7) 


2780.21(0.05) 


0.564(18) 


3.619(32) 


25.6 


f'24 


391.367(7) 


2555.15(0.05) 


0.525(18) 


4.925(34) 


23.9 


f 25 


679.150(30) 


1472.43(0.07) 


0.127(18) 


1.46(14) 


5.8 


ha 


977.465(37) 


1023.05(0.04) 


0.103(18) 


1.53(18) 


4.7 


hi 


2925.794(30) 


341.788(0.004) 


0.130(18) 


0.14(14) 


5.9 
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Figure 2b. Same as Fig. |2a]but for KIC 10001893, with the third panel 
showing residuals after removing 16 periodicities, and the bottom panel 
after 27 periodicities removed. The S/N = 4 threshold is at 0.088 mma. 



plitude spectrum of KIC 10001893, presented in the top panel of 
Fig. [2b] is also dominated by a number of peaks in the low fre- 
quency region. Two harmonics of the LC artefact (7 th and 8 th ) ap- 
pear in the high frequency range. Three remaining high-frequency 
peaks lie slightly above the detection threshold. Two of them are 
found in a "transition region" and one in the p-mode region. The 
latter, if real, might be a signature of hybridity. In the low frequency 
region we detected 24 peaks. The frequencies and amplitudes (of 
order of 1 mma) suggest that these peaks are associated with g- 
modes, so it is another V1093 Her star. The list of detected peaks is 
given in Table [4] 

The bottom panel of Fig.[2b]displays the final residual ampli- 
tude spectrum with all 27 detected peaks removed. The horizontal 
line represents a detection threshold of 4 times the mean noise level, 
which for this star is 0.088 mma. As is clearly seen, no peaks are 
left in the residuals at this level or above. This denotes that all peaks 
were resolved and no amplitude/phase variations on time scale of 
the run duration are present in this star. 

KIC 8302197 This star has the smallest number of peaks in 
the amplitude spectrum (Fig(2cJ in this sample. We found only 7 
peaks in the low frequency region, shown in the middle panel of 
FigHc] and one at high frequency with an amplitude above the de- 
tection threshold. All peaks are listed in Table [5] We can easily 
associate these seven peaks with gravity modes and identify this 
object as a V1093 Her-type star. The amplitude of the highest peaks 
is only 0.647 mma. 

The single peak at high frequency does not match any of the 
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Figure 2c. Amplitude spectra of KIC 8302197. The top panel shows the 
original amplitude spectrum up to the Nyquist frequency. The middle panel 
shows the relevant low-frequency domain, while the bottom panel shows 
the residual amplitude spectrum after removal of all 8 peaks. Short vertical 
dashed lines indicate the seven removed frequencies. The S/N = 4 threshold 
at 0. 1 64 mma is indicated by the dashed horizontal line in the bottom panel. 



known LC artefact frequencies, so we include it in our solution. 
Note that it was detected at S/N = 4 so its detection is marginal. 
Confirmation will require more data. As in the case of the previous 
objects, if its authenticity is confirmed this object will be another 
hybrid candidate. The amplitude spectrum also contains two strong 
LC artifacts (7 th and 8 th harmonics). 

The bottom panel of Fig[2c1 shows the final residual amplitude 
spectrum with all 8 detected peaks removed. This process did not 
reveal unresolved peaks with significant amplitudes. The time span 
of the data is sufficient to resolve all modes down to the detection 
threshold, which for this star is equal to 0.1 64 mma. 

KIC 11558725 The amplitude spectrum of this object, shown 
in Fig|2dl contains the largest number of detected peaks, listed in 
Table[6] A total of 41 peaks are revealed in the low frequency region 
between 50 and 400 /iHz. There are also quite a few in a "transition 
region" (9 peaks) as well as at high frequencies (3 peaks). Only two 
of them are found at high confidence. The S/N ratio for £43 and f.44 
is 12 and 17.4, respectively. It indicates that they are very signifi- 
cant and may be intrinsic to the star. We have to keep in mind, how- 
ever, that there is still chance that these peaks might be caused by 
a contaminating star; the contamination factor for KIC 1 1558725 is 
0.269 as currently listed in the Kepler Input Catalog. Unlike other 
objects presented in this paper, KIC 11558725 shows a significant 
number of peaks at higher frequencies. It is very likely that some 
(if not all) of them are real p-modes. This make the object best 
candidate for hybrid star. 

After removal of all 53 periodicities, the amplitude spectrum 
still contains residual signals, as can be seen in the bottom panel 
of Fig(2d] In particular, there are two strong peaks between 300 
and 350 /iHz. These two, along with other outstanding peaks above 
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Table 5. Results of the pre-whitening process for KIC 8302197. The phases are given at mean epoch BJD = 2455 169.508834. The numbers in parentheses are 
the errors of the last digits. 



ID 


Frequency [^tHz] 


Period[s] 


Amplitude [mma] 


phase [rad] 


S/N 


fi 


126.295(23) 


7918.0(1.4) 


0.354(33) 


5.3(0.8) 


8.6 


h 


145.588(12) 


6868.7(0.6) 


0.647(33) 


2.94(0.42) 


15.8 


f 3 


163.683(18) 


6109.4(0.7) 


0.443(33) 


5.5(0.6) 


10.8 


U 


178.335(46) 


5607.4(1.5) 


0.172(33) 


2.2(1.6) 


4.2 


fr, 


186.989(16) 


5347.91(0.46) 


0.497(33) 


1.0(0.5) 


12.1 


f„ 


283.243(33) 


3530.54(0.41) 


0.241(33) 


4.7(1.1) 


5.9 


h 


305.834(20) 


3269.75(0.21) 


0.409(33) 


0.1(0.7) 


10.0 


fs 


6002.352(48) 


166.601(0.001) 


0.165(33) 


4.8(1.6) 
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Figure 2d. Same as Fig.|2a]but for KIC 1 1558725. The third panel is the Fourier transform after 22 peaks are removed, and the bottom panel after 53 peaks 
removed. Short solid vertical lines mark the other 3 1 removed peaks. The S/N = 4 threshold is at 0.064 mma. 
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Figure 3. Best candidates for triplets in three stars. Vertical solid lines de- 
note the location of peaks which form the triplets. 

the horizontal dashed line at the detection threshold of 0.064 mma, 
were not removed, as they lie at frequencies within the frequency 
resolution limit of stronger peaks already removed from data. Ei- 
ther they are intrinsic to the star and caused by amplitude/phase 
changes, or result from beating between unresolved periodicities. 
With the longer data sets we expect to obtain, we hope to resolve 
this issue. We also found three harmonics of LC artifact (6* ,8* 
and 13 th ). 

KIC 10553698 This object also shows a very rich amplitude 
spectrum. As seen in Fig(2e] the spectrum is dominated by forest 
of peaks in the low frequency region characteristic of V1093 Her 
pulsators. We can easily identify at least 17 peaks with the highest 
amplitudes. Another 26 peaks (43 in total) with lower amplitudes 
were detected above the detection threshold of 0.076 mma and re- 
moved from data. We present our results in Table [7] 

Even after all these peaks have been pre-whitened, there is 
still residual signal present in the final residual amplitude spectrum 
(bottom panel in Fig[2et, The highest residual peak occurred very 
close to the dominant mode frequency. To sort out this problem we 
likely need more data to either resolve close frequencies or monitor 
amplitude/phase changes. 

Other than (/-modes, we found 4 signatures of LC artifacts 
(6 th , 7 th , 10 th and 14 th harmonics) and six unidentified peaks. 
Four peaks are located in the "transition region" and two are at 
high frequencies. All of them were detected with low significance 
and their authenticity has yet to be confirmed. 

In the bottom panel of Fig|2e]we can also notice a broad set of 
peaks at 360 /iHz. Although their amplitudes are above the detec- 



tion threshold and are not residual peaks of any formerly removed, 
we did not attempt to pre-whiten this region. This complex pattern 
is seen in many other compact stars observed with Kepler. It was 
seen in Q0 through Q4 data, and most frequently in the Q3 data re- 
lease. Although the cause of this feature is unknown and it is seen 
in a few stars, we assume that these peaks are caused by some in- 
strumental effects and were not included in our analysis. 

Comparing results for the two stars KIC 11558725 and 
KIC 10553698 we found that each star has a peak at nearly the 
same frequency: fsi in the former star and f±2 in the latter differ by 
0.4/iHz, which is comparable to the formal resolution of each run. 
We believe that these peaks are not intrinsic to the stars (either tar- 
get or contaminating ones) and likely to be caused by instrumental 
issues. Presumably f$2 in KIC 11558725 is a part of the complex 
structure and is also caused by the same instrumental effect. Since 
we cannot be absolutely sure about their origin yet we keep them 
in our solution. They do not influence peaks in the g-mode region. 



4 DISCUSSION 

With the available sets of frequencies for each star, we can now ex- 
amine them for the the asymptotic signature of (/-mode pulsations, 
the coherence of the observed pulsations, and the asteroseismic ef- 
fects of rotation and binarity. In this section, we identify the possi- 
ble signature of rotation (equally split multiplets), and discuss bina- 
rity (low-frequency photometric variation) and the possible hybrid 
behaviour of these stars. We also demonstrate the asymptotic sig- 
nature of g-mode pulsation, equal period spacings for consecutive 
overtones of the same degree I, but leave more detailed discussion 
of this property to a companion paper (Paper VIII). 

4.1 Rotationally split multiplets 

For nonradial pulsation in spherically stars, all modes with nonzero 
values of the azimuthal quantum number m have the same fre- 
quency as the m = mode. Azimuthal symmetry can be destroyed, 
lifting the rn degeneracy, by processes such as rotation or global 
magnetic fields. Since all stars rotate, we expect to see multiplets in 
the amplitude spectra, with the splitting within multiplets reflecting 
the rotation frequency. If all components of a multiplet are mani- 
fest, it will have 21+1 components. If the rotation is sufficiently 
slow (with a rotation period that is much longer than the oscilla- 
tion period) the central frequency (m = 0) will not be affected by 
rotation while all other components with non zero m values will be 
shifted by u n j^ m — v n> ifl + mtt(l — C„j) (Ledoux 1 95~0) in the 
simplest case of solid-body rotation at a frequency Q. The Ledoux 
parameter C n .i is an integral quantity that depends on the radial 
dependence of the oscillation eigenfunctions through the star. Iden- 
tification of multiplets in the amplitude spectrum can help to derive 
the rotation period, and is very useful for mode identification. 

Thu s far we detected multiplets in only a few V361 Hya class 
stars (e.g. lReed et al.l2 004; Baran et al. 2009). The very rich ampli- 
tude spectrum of V1093 Her class stars detected with Kepler hold 
the potential to expose multiplets in these stars. Without any obvi- 
ous evidence of these stars being in close binaries (with synchro- 
nized rotation), we do not have any prior expectation of what split- 
tings to look for (unlike the stars discussed in Paper V). Therefore, 
we constrained our search to rotation frequencies that range from 
the minimum that can be resolved with 30 days of data (on the 
slow rotation side) to rotation periods that would not produce ro- 
tational broadening of spectral lines greater than the limits placed 
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Table 6. Results of the pre-whitening process for KIC 11558725. The phases are given at mean epoch BJD = 2455 169.508302. The numbers in parentheses 
are the errors of the last digits. 



ID 


Frequency [/iHz] 


Period |s] 


Amplitude [mma] 


phase [rad] 


S/N 


h 


78.132(37) 


12798.9(6.1) 


O 087 
U.UO / 


1 3 1 
UJ 


1.1(6) 


5.5 


fa 


90.131(39) 


11095.0(4.8) 


O OR9 
U.UoZ 


UJ 


4.5(7) 


5.1 


f 3 


100.091(36) 


9990.9(3.6) 


O OQd. 


13) 


0.5(6) 


5.9 


£4 


101.002(33) 


9900.8(3.2) 


1 01 


1 3) 


4.6(6) 


6.5 


to 


104.524(24) 


9567.2(2.2) 


O 1 19 


1 ^ 1 

/ 


5.49(41) 


8.3 


fe 


116.106(24) 


8612.8(1.8) 


n 1 11 

U. 1 3 1 


1 1 \ 
13) 


0.31(42) 


8.2 


tr 


123.435(17) 


8101.4(1.1) 


O 1 87 


1-5 J 


4.04(29) 


11.7 


fa 


131.148(4) 


7624.97(0.23) 


O 8^9 


1 3 1 
lJj 


2.84(7) 


52.0 


f 9 


140.418(18) 


7121.6(0.9) 


O 1 80 


3) 


3.19(30) 


11.3 


fio 


145.076(30) 


6892.9(1.4) 


n 1 ns 

U. lUo 


13) 


5.48(50) 


6.8 


fn 


150.918(36) 


6626.1(1.6) 


O OQO 

u.uyu 


1-5 J 


3.7(6) 


5.6 


fl2 


152.612(19) 


6552.6(0.8) 


O I 71 




4.35(32) 


10.8 


fl3 


156.957(30) 


6371.2(1.2) 


O 1 07 


3) 


5.4(5) 


6.7 


fl4 


163.770(46) 


6106.1(1.7) 




13) 


6.2(8) 


4.3 


fir, 


170.378(24) 


5869.3(0.8) 


1 

U. iJJ 


1 3 1 

1-5 J 


1.58(41) 


8.4 


fl6 


176.204(10) 


5675.24(0.32) 


O 1A9 




3.47(16) 


21.4 


fl7 


178.229(37) 


5610.8(1.2) 


O 087 
\J.\)o 1 


1 3 1 
iJJ 


1.9(6) 


5.4 


fl8 


186.044(22) 


5375.1(0.7) 


142 


1 31 

IJJ 


3.96(38) 


8.9 


fl!) 


194.148(20) 


5150.7(0.5) 


161 


I 3 1 

1-5 J 


0.99(34) 


10.1 


fao 


203.747(5) 


4908.05(0.12) 


709 


1 31 

uj 


2.17(8) 


44.3 


I-21 


214.413(7) 


4663.90(0.15) 




i 3 1 

iJJ 


2.08(13) 


27.2 


f22 


226.192(8) 


4421.02(0.16) 






3.82(14) 


26.0 


f23 


226.996(20) 


4405.36(0.39) 


167 




5.88(34) 


10.4 


f24 


229.400(8) 


4359.20(0.15) 


425 


1 31 
1 J J 


2.29(13) 


26.5 


fas 


240.364(5) 


4160.36(0.09) 


O SQ'^ 

U.J7J 


1 31 
iJJ 


5.30(9) 


37.0 


f26 


251.902(26) 


3969.80(0.41) 


O 1 91 


IJJ 


3.46(44) 


7.7 


fa7 


261.621(45) 


3822.32(0.66) 


071 


I 3 1 

1-5 J 


0.4(8) 


4.4 


fas 


271.288(25) 


3686.12(0.34) 


129 


I 3 1 

1-5 J 


0.97(42) 


8.1 


fa9 


274.686(2) 


3640.52(0.03) 


1 19 1 
l.JZ 1 


I ^ 1 

13) 


5.450(41) 


82.5 


fid 


283.222(13) 


3530.80(0.16) 


0.248 


13) 


0.49(22) 


15.5 


fsi 


295.079(24) 


3388.92(0.28) 


0.137 


13) 


1.92(40) 


8.6 


f-i2 




3377 ?8fO 051 


0.930 


13) 


2 66(6) 


58 2 


f33 


305.200(26) 


3276.54(0.28) 


0.124 


13) 


0.05(44) 


7.8 


f 3 4 


307.604(10) 


3250.93(0.11) 


0.314 


13) 


1.34(18) 


19.6 


f 3 5 


311.336(10) 


3211.96(0.10) 


0.331 


13) 


5.77(17) 


20.7 


f«i 


320.842(7) 


3116.80(0.07) 


0.473 


13) 


4.82(12) 


29.6 


fi7 


322.493(5) 


3100.84(0.05) 


0.667 


13) 


4.00(8) 


41.7 


fix 


350.187(4) 


2855.62(0.03) 


0.749 


13) 


5.14(7) 


46.8 


f39 


364.732(42) 


2741.74(0.32) 


0.076 


13) 


2.0(7) 


4.8 


f4() 


382.609(8) 


2613.63(0.05) 


O 497 


1 ^ 1 

13) 


5.29(13) 


26.7 


f41 


390.899(16) 


2558.21(0.10) 


0.206 


13) 


6.16(26) 


12.9 


£43 


658.183(38) 


1519.33(0.09) 


0.084 


13) 


4.2(7) 


5.2 


f43 


710.810(17) 


1406.85(0.03) 


0.192 




1.47(28) 


12.0 


t'44 


712.361(12) 


1403.78(0.02) 


0.279 


13) 


5.26(20) 


17.4 


f46 


1039.500(50) 


962.00(0.05) 


0.068 


13) 


4.6(9) 


4.3 


f46 


1040.622(46) 


960.96(0.04) 


0.077 


13) 


1.3(8) 


4.8 


f47 


1041.727(38) 


959.94(0.04) 


0.089 


13) 


1.1(6) 


5.5 


f48 


1270.344(37) 


787.19(0.02) 


0.086 


13) 


3.1(6) 


5.4 


f49 


1370.516(44) 


729.65(0.02) 


0.072 


13) 


4.2(8) 


4.5 


fso 


1423.021(35) 


702.73(0.02) 


0.091 


13) 


4.4(6) 


5.7 


fsi 


3073.449(23) 


325.367(0.002) 


0.137 


13) 


0.45(40) 


8.6 


fsa 


3075.084(31) 


325.194(0.003) 


0.104 


13) 


3.3(5) 


6.5 


f 53 


4114.660(32) 


243.033(0.002) 


0.100 


13) 


3.2(5) 


6.2 



via spectroscopy. We assume that Ledoux parameter for p-modes 
(C n ,i =0.5, and C n ,2 =0.17 - see Paper V). These constraints limit 
us to splittings ranging from 0.825 to 13.5 /iHz, which correspond 
to rotation periods between 7 days and 10.3 hrs, respectively. We 
looked for triplets only, and we assumed that their symmetries can 



deviate within the formal frequency uncertainty, which is about 
0.02 /iHz. 

We do not find any potential candidates for triplets in 
KIC 7668647 or KIC 8302197. There are few possibilities in the 
other three stars. KIC 10001893 shows two triplet candidates. They 
are both slightly asymmetrical. One consists of fr, f$ and fg with 
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Figure 2e. Same as Fig.|2a]but for KIC 10553698. 17 peaks with the highest amplitudes were removed prior to producing the Fourier transform shown in the 
third panel, and 43 peaks were removed to produce the spectrum in the fourth panel. The S/N = 4 threshold is at 0.076 mma. 



splittings equal to 5.97 and 5.95 fiKz, while the other triplet is 
formed by fg, fio and fn with splittings equal to 7.30 and 7.33 /jHz. 
The fg peak appears in both triplets, which is clearly inconsistent 
with a simple rotational splitting inteipretation. Since the amplitude 
of the modes in the 5.96 /^Hz-spaced triplet are significantly larger, 
we have more confidence that it is a true multiplet. 

In KIC 11558725 we found two other candidates. These are: 
fis, fi7 and fis, having splittings of 7.85 and 7.82 /jHz, and f45, f46 
and f47 with splittings equal to 1.12 and 1.11 /jHz. These clearly 
are not consistent with a uniform rotation period within the star, and 
even models with diff erential rotation (with d epth) rarely show this 
extreme difference dKawaler & Hostlerll2005l) . For KIC 10553698 
we find a candidate triplet formed by fg, fin and fix, with splittings 
of 6.17 and 6.19 /jHz. The amplitude spectra surrounding these 
triplets are presented in Fig. [3] 

Although we found a few candidates for rotationally split 
triplets, it is hard to explain why only very few modes within a 
given star are rotationally split. But there are similarities with other 



sdBV stars. For example. iBaran et"aL 1 d2009h found only two mul- 
tiplets among more than fifty periodicities in Balloon 090100001, 
and those involved the largest-amplitude modes. Assuming we are 
seeing real triplets here, a rough estimate of the average rotation 
period for those three stars is about 20 hrs. 



4.2 Binarity 

iMaxted et ail J200lh and lMorales-Rueda et all d2006t) . among oth- 
ers, argue that approximately half of all sdB stars are members 
of binary systems. A signature of binarity can be found in pho- 
tometric observations (via eclipses, modulation from reflection 
off of a low-mass companion, and/or through ellipsoidal varia- 
tions) or through radial velocity variations measured via spec- 
troscopy. ^larry_bmary_s^B_s^arsdrtected in p hotometry, e.g. 
HW Vir dWood. Zhang & Robinsorl 19931) , NY Vir dKilkennv et al] 
Il998h having beau tiful eclipses features or KPD 1930+2752 
i Billeres et all [2000) showing ellipsoidal variation. Recently an- 
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Table 7. Results of the pre-whitening process for KIC 10553698. The phases are given at mean epoch BJD = 2455200.902734. The numbers in parentheses 
are the errors of the last digits. 



ID 


Frequency [/iHz] 


Period |s] 


Amplitude [mma] 


phase [rad] 


S/N 


fl 


104.291(35) 


9588.6(3.2) 


0.094 




3.76(17) 


4.9 


h 


106.562(27) 


9384.2(2.4) 


123 


1 f\\ 

1 OJ 


3.29(13) 


6.5 


h 


109.701(35) 


9115.7(2.9) 


0.095 


1 (Si 

1 Of 


3.60(17) 


5.0 


U 


116.395(23) 


8591.4(1.7) 




1 (Si 
LOJ 


1.15(11) 


7.7 


to 


120.037(32) 


8330.8(2.2) 


104 


LOJ 


1.79(15) 


5.5 


fe 


128.775(39) 


7765.5(2.4) 


084 


1 0) 


1.00(19) 


4.4 


t'r 


137.569(7) 


7269.08(0.37) 


475 


1 (VI 
LOJ 


3.869(34) 


25.0 


f 8 


140.513(4) 


7116.78(0.20) 


842 


1 (Si 

LOJ 


0.973(19) 


44.3 


f<> 


148.313(8) 


6742.50(0.36) 


434 


1 (Si 
LOJ 


1.964(37) 


22.8 


fio 


154.486(27) 


6473.1(1.1) 


121 


1 (Si 


4.91(13) 


6.4 


fn 


160.672(16) 


6223.9(0.6) 


201 


1 0) 


4.08(8) 


10.6 


fl2 


167.904(21) 


5955.8(0.7) 


154 


1 (SI 
IOJ 


6.19(10) 


8.1 


fl3 


175.342(12) 


5703.14(0.39) 


9RS 


IOJ 


5.72(6) 


15.0 


fl4 


185.728(34) 


5384.2(1.0) 




1 (Si 
IOJ 


3.40(17) 


5.1 


fir, 


192.202(34) 


5202.9(0.9) 


098 


1 (Si 
IOJ 


0.44(16) 


5.1 


fl6 


201.963(2) 


4951.40(0.05) 


1 473 


1 (Si 

LOJ 


0.780(11) 


77.5 


fl7 


208.060(8) 


4806.31(0.18) 
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other binary sdB stars have been found using Kepler data 
dKawaler et ai]|2010bl ; l0stensen et a"i]|2010cl) . The first describes 
two pulsating sdB stars that are reflection-effect binaries; the sec- 
ond is a pulsator in an eclipsing binary. 

Analysis of the five stars presented here does not reveal any 
evidence for either the reflection effect or ellipsoidal variation or 
eclipses. We note that the de-trending described in Sect. 2 might 
mask some of that evidence, so we also analysed the data prior 
to de-trending. While detailed radial velocity time series have not 
been performed, there is no evidence for composite spectra in these 
targets. Therefore, we have no evidence of binarity for these stars. 
We note that with a longer time baseline, timing variations of the 
principal pulsation modes can be us ed to reveal a com panion via 
reflex orbital motion, as was done bv lSilvottiet"al1 l l2007l) to reveal 
a planetary companion to V391 Peg. 



4.3 Do these stars show p-mode pulsations? 



As mentioned in Sect. 1, the two classes of pulsating sdB stars 
overlap in the H-R diagram. On the hotter end, we find p- 
mode pulsators (the V361 Hya class). Typically, these are hotter 
than T c ff ~30kK. On the cool side we have g-mode pulsators 
(V1093Her) with typical T cft of about 25 kK. In the overlap re- 
gion, stars can show both g-mode and p-mode pulsations. Since 
these modes sample different regions of the stars, they are particu- 
larly potent targets for asteroseismic investigation. Prior to the Ke- 
pler discoveries, these stars were relatively rare, with only 5 ex- 
amples found in ground-based studies. In all of those cases, the 
stronger pulsations are in the p-mode frequency range. The lone Ke- 
pler p-mode pulsator had a single low-amplitude oscillation at low 
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frequencies (Paper II), suggesting that it behaves like the already- 
known hybrid pulsators. 

For the Kepler sdBV stars from the first part of the survey, we 
found that many of the <?-mode pulsators showed suggestive evi- 
dence for shorter period pulsations in the p— mode regime (Paper 
III) at lower amplitude. This pattern is not seen in any of the sdBV 
pulsators observed using ground-based data. This is likely a signal- 
to-noise issue. The g-modes are low amplitude to begin with, and at 
problematic frequencies for adequate ground-based observing; the 
possible p-modes we see are at even lower amplitude. 

Based on their effective temperatures, almost all the potential 
Kepler hybrid candida tes (five stars in this paper together with two 
stars from paper by iReed et al]|2010l) have a typical T e ff close to 
27 kK. The only exception is the slightly-cooler KIC 5807616 with 
only one barely detected peak at high frequencies. Their location in 
log g and T e g diagram is right between hybrid stars dominated by 
p-modes and V1093 Her class stars. This would be consistent with 
the idea that there is a transition temperature between p-mode pul- 
sations at the hot end and g-modes at the cool end; the p-mode dom- 
inant stars are in general hotter than the Kepler (/-mode dominant 
hybrids. If the single high-frequency periodicities are confirmed, 
this could explain their presence. 

The issue of isolated p-modes is addr essed in more detail in 
one of these stars in the modeling work by Charpinet et al.l ( 20 1 1 ) 
of the star KIC 2697388. Their theoretical models suggest that one 
peak at high frequency is likely intrinsic to the star, despite its low 
amplitude. Longer coverage should provide a definitive test of the 
reality of this mode. 



4.4 Amplitude and phase stability 

While pre-whitening we assumed that each periodicity maintained 
a fixed frequency, amplitude and phase. This does not have to be 
the case, however, since the pulsation can experience amplitude 
and phase modulation through nonlinear interactions, or growth or 
decay of mode energy over the course of the run. Changes in ampli- 
tude and/or phase can produce a broadening of the peaks in the am- 
plitude spectrum, or produce multiple closely-spaced peaks. Mani- 
festation as multiple peaks could mimic the appearance of multiple 
(fully coherent) pulsation modes, so resolving the amplitude spec- 
trum requires care to distinguish between these possibilities. 

For some stars, pre-whitening leaves behind no signals above 
our noise limits in the residual amplitude spectra. When no resid- 
ual signal at the frequency of a removed peak remains, we can be 
sure that the peak was resolved and the frequency, amplitude, and 
phase are constant over the course of data run. Looking back at the 
residual spectra obtained in our analysis we can see that no residual 
signals remain after pre-whitening the data from KIC 8302197 and 
KIC 10001893. Removal of the tabulated peaks leaves no signifi- 
cant signal behind. 

However, for three stars, residual peaks remain after pre- 
whitening the dominant periodicities. Their residual spectra con- 
tain additional signal at the frequencies of formerly removed peaks. 
These may be closely spaced coherent modes, or the effects of am- 
plitude and phase modulation. For KIC 7668647 this extra signal 
remains in four places with relatively small amplitudes. For two 
other stars, KIC 1158725 and KIC 10553698, significant peaks do 
remain after prewhitening. In the former, there are two strong peaks 
slightly above 300 pHz while there is one strong residual peak in 
the latter. In both cases there are other peaks with smaller signifi- 
cance but still above the detection threshold. 
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Figure 4b. Same as Fig.ga]but for KIC 1 1558725. 



There is a correlation between the amplitude and phase varia- 
tions that is characteristic of beating between closely-spaced coher- 
ent oscillations (i.e. the phase jumps by 180° during the amplitude 
minimum in a beat between two equal-amplitude periodicities). If 
the amplitude and/or phase vary irregularly and with no specific 
correlations, then that might suggest that the modes are experienc- 
ing stochastic changes. Since we have continuous data that span 
about 30 days we can use shorter subsets of the data to look for 
amplitude and/or phase variations, as was done in Paper II. 

To examine the phase and amplitude stability at a given fre- 
quency, we first prewhitened by all of the other significant peri- 
odicities found for that star. The prewhitened data were then di- 
vided into shorter, equal-length, contiguous runs. For each of these 
shorter runs, we fixed the frequency at the best value for the mode 
in question, and fitted the amplitude and phase. Choosing the length 
of subsets requires balance between time resolution and amplitude 
and phase uncertainties. We performed this analysis with different 
subsets, and found that 2 day samples provided the best compro- 
mise. We used non-linear least-squares to fit one sine wave to each 
individual subset keeping frequency fixed for each peak to deter- 
mine the pulsation amplitude and phase (using a fixed but arbitrar- 
ily zero point). We repeated this fitting only for four highest resid- 
ual peaks for each of those three stars. The result of this analysis is 
presented in Fig |4al4cl 

Fig. l4al shows that the amplitude of this star changes at a sig- 
nificant level for all four tested frequencies. In all of these fre- 
quencies, there is no clear correlation between the amplitude and 
phase modulations. For example, there are no abrupt phase jumps 
in any of the modes tested; though phase variations occur, they are 
all much less than one cycle, and they are not correlated with a 
minimum in pulsation amplitude. Thus the variations appear to be 
caused by random variation in the parameters of this mode. The 
variations seen in three of the four modes shine in Fig. |4b|appear 
to be stochastic as well, though in fzi there appears to be a jump in 



phase around day 167, when that mode is at minimum amplitude. 
The period of this mode is 3250 seconds, and the phase jump is 
about 1600 seconds - suggesting that there is beating present near 
this mode. 

KIC 10553698 shows much slower, smoother variations in 
amplitude and phase in the four modes tested as compared to the 
other two stars. In most cases, there is no correlation between the 
phase and amplitude variations for a given mode (i.e. no large phase 
changes near amplitude minima). But as can be seen in Fig. [4c] ft 
shows an unusual phase variation near minimum amplitude around 
day 200-205. The period of this mode is in excess of 7000 seconds, 
and the phase increases in that interval by 1500 seconds and then 
returns to the slower trend. It is conceivable that the phase jumped 
by half a cycle (3500 seconds) and that the (O-C) curve should lie 
up near 3500 seconds in this plot if beating is present. But our best 
guess is that the phase did not experience that large a change. It is 
notable that the slow variation in the amplitude of fie is reflected in 
the large residual peak near 200 /iHz in the bottom panel of Fig[2el 
Clearly, a longer observing run on this star should aid in determin- 
ing the true nature of these periodicities. 



4.5 Asymptotic (/-mode period spacings 

Theo retical consideration of the VI 093 Her class dFontaine et al.l 
l2003h showed that the pulsation modes detected in these stars cor- 
respond to high order (/-modes. As such, we could invoke asymp- 
totic limit for n > ! in which consecutive overtones should be 
evenly spaced in p eriod. This behavior has be en observed in pulsat- 
ing white dwarfs dKawaler"& Bradlev 1994). Because of the low- 
amplitude signal that is characteristic of these stars, and the logisti- 
cal challenges of observing these frequencies from a single ground- 
based observing site, this asymptotic behaviour had not been seen 
in VI 093 Her stars until space-based observations became avail- 
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able. The first set of Kepler VI 093 Her stars were found to show 
clear evidence of equal period spacings, with values of approxi- 
mately 145 and 250 s reported in Paper III and in a later work (Pa- 
per VIII). 

The five new detected VI 093 Her class also show evidence 
of evenly-spaced periods. For all stars we found average spacing 
in 250-270 s range. According to the Paper III, all of the peaks 
matching this spacing are likely to be I = 1 modes. There is also ev- 
idence of a second series of peaks with spacings in the 140-150 s 
range, consistent with them being I = 2 modes. This simple relation 
between consecutive overtones provides a determination of the de- 
gree I of the modes, which will be crucial for interior modeling 
efforts. Detailed analysis of period spacing on five stars presented 
here and stars already published in previous papers of our series is 
a subject of Paper VIII. 



5 SUMMARY 

We found five new long-period pulsating sdB (V1093 Her) stars 
using the Kepler photometer. The oscillations are mostly in the 
100^-00 /.tHz region, although we also found some single peaks 
at higher frequencies. The nature of those higher-frequency peaks, 
which are at low amplitude, remain uncertain, and we cannot rule 
out contamination from variable contaminating sources, or un- 
known spacecraft artefacts. On the other hand, assuming that peaks 
at high frequencies are indeed p-modes, the stars match the previ- 
ously known T c fr - pulsation type correlation. We believe that more 
data should confirm or disprove our hypothesis and definitively an- 
swer the question about the origin of the peaks. The amplitudes of 
detected peaks in the low frequency region are about 1% or less 
which is in similar to other known V1093 Her stars. 

As we have not found any signatures of binarity we claim that 
these are single stars. The pulsation spectra show only limited evi- 



dence of rotational splitting, but these occur with low-amplitude pe- 
riodicities and therefore require continued observation to improve 
the signal-to-noise and look for other similar spacings in the same 
stars. While many of the observed oscillations are coherent, sev- 
eral modes show the signature of possible stochastic amplitude and 
phase modulation. 

This paper completes our initial analysis of the survey phase 
data on compact stars by Kepler. In total we found 12 new mem- 
bers of the VI 093 Her class (5 in this paper, 2 in Paper V, and 
5 in Paper III), one V36 1Her star (Paper II) and one hybrid 
star l l0stensen et alj|2010ch . Among the former two groups, many 
showing additional peaks in the opposite side of the amplitude 
spectrum. 

Our results demonstrate the potential of the Kepler spacecraft 
for asteroseismology of sdB stars. Continuous coverage and the un- 
precedented quality of the data allowed us to detect, at high S/N, 
stellar pulsations in stars which would not be easily studied from 
the ground. The rich amplitude spectra of V 1093 Her class stars en- 
ables us to use the relation for period spacings to assign most of the 
detected peaks with n (relative) and I values (see Paper VIII). With 
those determinations, along with the frequencies themselves, tight 
constraints can be placed on models of these stars, and therefore 
will expand our understanding of the formation, structure, and evo- 
lution of these highly evolved stars. We anticipate that with much 
longer runs of Kepler data we can detect even more peaks with 
lower amplitudes making period spacing tool even more efficient. 
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